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Baryon history and cosmic star formation in non-Gaussian 
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ABSTRACT 

We present the first numerical, N-body, hydrodynamical, chemical simulations of cosmic 
structure formation in the framework of non-Gaussian models. We study the impact of pri- 
mordial non-Gaussianities on early chemistry (e~, H, H + , H~, He, He + , He ++ , H2, H^, D, 
D + , HD, HeH + ), molecular and atomic gas cooling, star formation, metal (C, O, Si, Fe, Mg, 
S) enrichment, population III (popIII) and population II-I (popll) transition, and on the evolu- 
tion of "visible" objects. 

We find that non-Gaussianities can have some consequences on baryonic structure formation 
at very early epochs, but the subsequent evolution at later times washes out any difference 
among the various models. When assuming reasonable values for primordial non-Gaussian 
perturbations, it turns out that they are responsible for: (i) altering early molecular fractions 
in the cold, dense gas phase of ~ 10 per cent; (ii) inducing small temperature fluctuations 
of < 10 per cent during the cosmic evolution of primordial objects; (iii) influencing the on- 
set of the first star formation events, at z > 15, and of the poplll/popll transition of up to 
some 10 7 yr; (iv) determining variations of < 10 per cent in the gas cloud and stellar mass 
distributions after the formation of the first structures; (v) causing only mild variations in the 
chemical history of the Universe. We stress, though, that purely non-Gaussian effects might 
be difficult to address, since they are strictly twisted with additional physical phenomena (e.g. 
primordial gas bulk flows, unknown primordial popIII stellar mass function, etc.) that have 
similar or stronger impact on the behaviour of the baryons. 
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1 INTRODUCTION 

Cosmic structures are supposed to rise from primordial matter fluc- 
tuations generated at very early times, during the inflation era. 
These fluctuations represent the seeds that grew over the cosmic 
time till the formati on of the presently observed Universe (e.g. 
iKomatsu et alj|201~ih . In the standard scenario, a Gaussian distri- 
bution of such perturbations is assumed, as a consequence of the 
central limit theorem. 

However, cosmic microwave background (CMB) temp erature fluc- 
tuations, at high-order perturbation theory (see e.g. IZaldarriagal 
200d ; lBartolo et alJ|2004l:ISeerv & Lidsevl|2005l: lAssadullahi etaij 
20071: ICoorav et al.Il2008l: iBeltranl 120081: iLiguori & Riottol 120081 



ulations of large-scale structures have shown how the abun- 
dance of rare-peak dark-matter haloes and thei r clustering prop- 



erties and distribution migh t be affected (e. g. Grinstein & Wise! 
1 19861 : |Kovamaet alJ 1 19991; iRobinson et aljl200d: IGrossi et alj 

iGrossi etalTT2009: 



20071: kang et alj 12001 iDalal et all 120081: 
Pesiacaues et alJl2009l ; lLam et alfeOldlMaggiore & Riottdl201 
I Pillepich et alj |201C|: iRoncarelli et alj l20ld : IPace et alj 1201 



2011 



Khatri & Wandeld |2009l: |Lam & Shetlj 120091: iBartolo et al 



2010; 



Yadav & Wandeltl |2010| : iGad l2010l : [Pitrou et al l l2010h . could 
have deviations from the purely Gaussian shape and leave 
room for non-Gaussian models. Additional studies of the ex- 
cursion set formalism and analyses of N-body numerical sim- 
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Wagner et alj \20 lot bo Verde & Smith! 1201 li iMaturi et alj 
lYokovamaet all201lh . ~ 
Observationally, the CMB is the best way to probe non- 
Gaussianities and different experiments (COBE, BOOME RanG , 
WMAP, PLANCK) can be used (e.g.lKomatsu et aT]|200ll200j ; 
Gaztanaga & Wagdl2003l : [Spergel et aI1l2007l : iHikage et alj|2008l: 
Yadav & WandelJ T2OO8F lAfehordi & Tollevl I2OO8I: Isiosar etal] 
20081: IKomatsu et alj 120091: iNatoli et alj 120091: iKomatsij |20K 



Hou et alj l20ld : iRath et alj 120091 : iRaeth et alj l20ldl . Recent 
determinations of the cosmological parameters by the 7-year 
Wilkinson Microwa ve Anysotropy Probe (WMAP) satellite (e.g. 
IKomatsu et alJt20Tlh still find primordial non-Gaussianities. 
Despite the huge amount of work based on N-body dark-matter 
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only simulations, detailed chemical and hydrodynamical investiga- 
tions of the effects on the visible matter are still missing, or poorly 
under stood. According to previous numerical studies dViel et al.l 
2009), the Lyman-a flux bispectrum in the high-redshift intergalac- 
tic medium is expected to present deviations from the Gaussian 
case up to ~ 10 pe r cent at z ~ 4. Semianalytic estimates (e.g. 
ICrociani et afl |2009) seem to show that the ionized fraction of the 
intergalactic medium (IGM) can grow by a factor of a few, up to 
5, with respect to the corresponding Gaussian model. The increase 
of the filling factor has a small impact on the reionization optical 
depth and is of the order of ~ 10 per cent if a scale-dependent 
non-Gaussianity is assumed. Since non-Gaussianities are expected 
to influence the primordial matter density fluctuations, it is impor- 
tant to address their consequences on the history of the baryons, 
from the early "dark ages", when stars were not formed, yet, to the 
following "louminous ages", after the birth of the first objects. Ad- 
ditionally, it is relevant to understand the role they play in the transi- 
tion from the primordial population III (popIII) stars to present-day 
population II-I (popll) stars. These are indeed among the goals of 
the present work. 

A standard way to account for the existence of non-Gaussianities 
in the primordial fluctuation field is to regard them as a pertur- 
bation on a Gaussian background. Strongly non-Gaussian models 
have been progressively dismissed by CMB observations and it is 
now generally agreed that primordial non-Gaussianities, if present, 
must have represented a moderate deviation from the underlying 
Gaussian statistics. The most commonly adopted parametrisation 
of non-Gaussianities considers them as second-or der perturbations 
of the Bardeen gauge-invariant potential (see e.g.lSalopek & Bone 
19901 : IKomatsu & Spergelll200ll : IVcrddbOlOb iDesiacques & Seliak 



2010) 



* = *L + /NL [<&l- < *L >] 



(1) 



where $l is the linear Gaussian part, and the dimensionless cou- 
pling constant, /nl, controls the importance of the non-Gaussian 
contributiorQ. Since the final value of the potential $ depends 
on the local value of the Gaussian field $l, non-Gaussianities of 
the the form ([]} are labelled as "local". Local models are physi- 
cally relevant for scenarios where non-linearities develop outside 
the horizon, but they are by no means exhaustive of all the forms 
non-Gaussianities can assume according to the several mechanisms 
that generate them (for further details se e e.g iBabich et al.ll2004l ; 
Bartolo et al. 20041: ILo Verde et"al]|2008l and references therein). 
The exact value of the /nl parameter is continuously better con- 
strained by CMB observations (see discussion in Sect. |4). More 
extended analyses could in principle be done by introducing addi- 
tional parameters, like the cubic-order one, <7nl, but, given their 
large uncertainties (e.g. —7.4 x 10 s <jnl< 8.2 x 10 5 , accord- 
ing to lSmidtet afll2010l) . we will consider only second-order cor- 



The Bardeen gauge-invariant potential, reduces to the usual New- 
tonian peculiar gravitational p otential, up to a minus sign, on sub-Hubble 
scales. In the literature (e.g. Verde 2010), there are two, sometimes mis- 
leading, notations: the large-scale structure (LSS) convention and the CMB 
convention. In the former (LSS), <E> is linearly extrapolated at z = 0; in 
the latter (CMB), $ denotes the primordial potential, at z — > +oo. The net 
difference is related to the ratio of the growth factor, g(z), at the two times 



and accounts for a correction of roughly unity: i.e. /n 



9(z 



-oo)/g(z = 0) / NL CMB , which corresponds to ~ 1.28/ NL CMB , for 
a standard ACDM model, a nd to ~ 1.33/nl CM B , for the latest obser- 
vational data from WMAP ( Komatsu et al, 2011). Usually, observational 
values are reported according to the CMB convention. 



rections. In addition , since latest determinations suggest /nlS", 
i Spergel et al.l 120071: Isiosar et al.l [20081: lAfshordi & Tollev) 2008; 
Yadav & Wandeltl l2008t ICurto et al.l 120091: IKomatsu et alj I2OO9I; 
Smith et al.ll2009l;lRudiord et alj|2009t lKomatsull2010l : iHou etal] 
201fJ : lRudiord et al.ll2010l: ICabella et alj|201fj|) . we will only con- 



sider positive values for /nl. 

In the present work we will focus on the effects non-Gaussianities 
have on the baryon history in universes with different initial con- 
ditions, by studying how gas evolution, cooling and condensa- 
tion could lead to different star formation epochs for different ini- 
tial scenarios, and by exploring how the various stellar population 
regimes and the statistical properties of gas and stars are affected in 
non-Gaussian models. 

We will present the simulations performed in Sect. [2] together 
with the main results, in Sect. [3] about the first cosmic struc- 
tures (Sect. I3.lt . star formation (Sect. 13.2) . chemistry evolu- 
tion (Sect. 13.3b and primord i al streaming mo tions (see also e.g. 
Tseliakhovich & Hirata 2010; Maio et aU201ll) . and baryon statis- 
tics (Sect. [3~l4t . Then, we will discuss and conclude in Sect.|4] 



2 SIMULATIONS 

We use a modi fied version o f the parallel tree/SPH numerical 
code P-Gadget2 ( Springel 2005). Beyond gra vity and hydrodynam- 
ics, it includes radi ative gas cooling a t high dSutherland & Dopita 
1 19931) and low (Ma io et al. I l2007l) temperatures, both from 
molecules and fine-structure metal transitions, multiphase model 
( Springel & Hernquist 2003) f or star formation (inspired on the 
works by iKatz & GunrJ Il99ll ; ICer] [T99l : ICen & Ostriken 1 19921; 
Katj Il992l: iKatz et all 1 19921 , Il99rj ), UV background radiation 
( Haardt & Madau|1 996l), a nd wind feedback dSpringel & Hernquistl 
2003l : lAguirre et al.ll200ll) . 



The main agents of gas cooling, in hot environments, are atomic 
resonant transitions of H and He excitations that can rapidly 
bring the temperatures down to ~ 10 4 K. In metal-rich environ- 
men ts, atomic fine-st ructure transitions can further lower them 
(see iMaio et alj |2007). mo stly at high redshift, since the forma- 
tion of a UV background dHaardt & Madaulll996l; lHaehnelt etail 
l200ll : iGilmore et alj [2009) at lower redshifts competes by heat- 
ing the IGM. In dense environments, stochastic star formation is 
assumed, cold gas is gradually converted into stars (that repro- 
duce the Kennicut-Schmidt low), outflows take place, and feed- 
back effects inject entropy into the surrounding environment, lead- 
ing to a self-regulated star forming regime (for further detail see 
ISpringel & Hernquistll2003l) . 

In the code, we also include the relevant chemical network to self- 
consistently follow the cosmic evolution of e — , H, H , H~, He, 
He+ He ++ , H 2 , H+, D, D+ HP HeH+ (see e.g. lYoshida etail 
l2003l:lMaio et al J 20071. 1 20091. 1 20 1 Ol and references therein), metal 
(C, O, Si, Fe, Mg, S) pollution from popIII and/or popll stellar 
generation s, ruled by a cr i tical metallicity threshol d of Z cr u = 
10~ 4 Z Q dTornatore etal] 120071; IMaio etal] |2010|) a nd leading 
gravi tational enrichment into the surrounding medium (IMaio et al.l 
l2010h . 

We perform several runs (see Tab.Q] in order to check the different 
assumptions on /nl, at large and small scales. 
The initial conditions are generated with a modified version of the 
N-GenIC code that introduces non-Gaussian features of the lo- 
cal type in the initi al density field according to eq. ifl) (for fur- 
ther details see e.g. iBartolo et al]|2005l : iGrossi et al.l 12007. 2009; 
IVielet alj [2009). We assume a concordance ACDM model with 
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Table 1. Initial parameters for the different runs. 



Runs 


Box side 


Particle mass [Mq//i] for 


Softening 


./nl 


Gas bulk 


PopIII IMF 




[Mpc/ft] 


gas (dark matter) 


[kpc//i] 




velocity [km/s] 


range [M ] 


Run05.0 


0.5 


42.35 


(275.28) 


0.04 








[100. 500] 


Run()5.10 


0.5 


42.35 


(275.28) 


0.04 


10 





[100, 500] 


Run05.50 


0.5 


42.35 


(275.28) 


0.04 


50 





[100, 500] 


Run05.100 


0.5 


42.35 


(275.28) 


0.04 


100 





[100, 500] 


Run()5.1000 


0.5 


42.35 


(275.28) 


0.04 


1000 





[100, 500] 


Run 100.0 


100 


3.39 x 10 s 


(2.20 x 10 9 ) 


7.8 








[100, 500] 


Run 100. 10 


100 


3.39 x 10 8 


(2.20 x 10 9 ) 


7.8 


to 





[100, 500] 


Runl00.50 


100 


3.39 x 10 s 


(2.20 x 10 9 ) 


7.8 


50 





[100, 500] 


Run 100. 100 


100 


3.39 x 10 8 


(2.20 x 10 9 ) 


7.8 


100 





[100, 500] 


r\ 1 f\r\ 1 f\r\f\ 

Run 100. 1000 


100 


3.39 x 10 8 


(2.20 x 10 9 ) 


7.8 


1000 





[100, 500] 


RunlOO.O.SL 


100 


3.39 x 10 s 


(2.20 x 10 9 ) 


7.8 








[0.1, 100] 


RunlOO.10.SL 


100 


3.39 x 10 8 


(2.20 x 10 9 ) 


7.8 


to 





[0.1, 100] 


RunlOO.50.SL 


100 


3.39 x 10 s 


(2.20 x 10 9 ) 


7.8 


50 





[0.1, 100] 


RunlOO.100.SL 


100 


3.39 x 10 8 


(2.20 x 10 9 ) 


7.8 


100 





[0.1, 100] 


RunlOO.1000.SL 


100 


3.39 x 10 8 


(2.20 x 10 9 ) 


7.8 


1000 





[0.1, 100] 


Run05.0.30 


0.5 


42.35 


(275.28) 


0.04 





30 


[100, 500] 


Run05.0.60 


0.5 


42.35 


(275.28) 


0.04 





60 


[100, 500] 


Run05.0.90 


0.5 


42.35 


(275.28) 


0.04 





90 


[100, 500] 


Run05. 100.30 


0.5 


42.35 


(275.28) 


0.04 


100 


30 


[100, 500] 


Run05. 100.60 


0.5 


42.35 


(275.28) 


0.04 


100 


60 


[100, 500] 


Run05. 100.90 


0.5 


42.35 


(275.28) 


0.04 


100 


90 


[100, 500] 



matter density parameter flo.m = 0.3, cosmological density pa- 
rameter firj.A = 0.7, baryon density parameter f2rj,b = 0.04, ex- 
pansion rate at the present of Ho = 70km/s/Mpc, power spec- 
trum normalization via mass variance within 8 Mpc//i radius sphere 
erg = 0.9, and spectral index n = 1. We study the non-Gaussianity 
parameter space identified by /nl= 0, 10, 50, 100, 1000, in cos- 
mological periodic boxes of 0.5 Mpc//i and 100 Mpc//7 a side, 
sampled, at redshift z — 100, with 2 x 320 particles with gas 
mass resolution of ~ 40M Q //i and ~ 3 x 10 8 M Q //i, and cor- 
responding maximum physical softening lengths of 0.04 kpc/A and 
7.8 kpc//j, respectively. This will allow us to see the impacts of 
the different models down to scales of the order of ~ 0.1 kpc. 
We note that the small sizfl of the 0.5 Mpc/h boxes prevents 
the simulations from being run to very low redshifts, in order for 
the fundamental mode to remain linear; the large sim ulations sat - 
isfy instead this requirement (see e.g. Sect. 2.3 in lMaioll2009l) 
and are thus run until z = 0. The simulations will be labelled 
as Run05.0, Run05.10, Run05.50,Run05.100, Run05.1000, and 
RunlOO.O, RunlOO.lO, Runl00.50,Runl00.100, Runl00.1000, re- 
spectively. 

Throughout this work, we will assume a Salpeter IMF for popll 
star forming regions (where Z ^ Z cr i t ), and a top-heavy IMF, 
with mass range [lOOM©, 500M©], for popIII star forming re- 
gions (where Z < Z cr u). The latter assumption is matter of 
debate, since there are also theoretical and numerical evidences 



2 Results on biased high-density (zoomed) initial conditions, for which 
the onset of star formation is expected to ta ke place at redshi fts as high 
as ~ 40 — 50 heve been presented by e.g. iMaio et at] fc009l) . While 111 
such set-ups the timing of the popIII formation might b e more accurate and 
independent from Lym an-Werner radiation feedback (John son et alj|200^ ; 
iTrenti & Stiavelli 2009), the general physical trends are supposed to be very 
similar. 



for the existence of low-mass popIII stars (e. g. Yoshidal 20061: 
l Yoshida et al.l200llCampbell & Lattanziol2008MSuda & Fuiimotol 
l2010l : IStacv et alj2010h . If this is the case, that might mean that the 
conclusions of several previous works predicting massive primor- 
dial stars were highly affected by numerics and resolution. As the 
primordial popIII IMF can impact the cosmological pollution his- 
tory (e.g. lMaio et ai] |2010h , in Sect. 13.21 we we will also consider, 
for the 100 Mpc//? side boxes, the opposite extreme of a Salpeter- 
like popIII distribu tion, with ap propriate yields, for primordial stars 
(see discussion in IMaio et al.l I2OI0L and references therein), and 
mass range [0.1M©, lOOM©]: the corresponding runs will be la- 
belled by adding the SL suffix. 

In order to see how relevant the consequences of non-Gaussian 
perturbations are, in Sect. 13.21 we will also compare them 
with second-order effects in the linearized perturbation theory 
dTseliakhovich & Hiratal l2010h . that expect supersonic gas bulk 
flows at early times to suppress structure growth, and delay star 
formation on scales smaller than a few Mpc dMaio et al .1201 ll) . We 
will focus on the /nl=0 and / NL =100 cases, assuming primordial 
bulk shifts of 30, 60, 90 km/s. The six corresponding simulations 
will be identified by Run05.0.30, Run05.0.60, Run05.0.90 (for the 
/nl=0 set), and by Run05. 100.30, Run05. 100.60, Run05. 100.90 

(for the / NL =100 set), respectively. 

A friend-of-friend (FoF) algorithm (Springel et al. 200 lj), with co- 
moving linking lenght of 20 per cent the mean inter-particle sepa- 
ration, is applied at postprocessing time, considering all the types 
of particles, to find the formed cosmic objects, with their dark, 
gaseous, and stellar components. 

We summarize the properties of all the different runs in Table[T] 
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3 RESULTS 

In the following, we will show the main repercussions primordial 
non-Gaussianities have on the cosmic baryon history in the Uni- 
verse. We will focus on cosmic structure evolution (Sect. 13. IK 
star formation (Sect. 13.2) . chemical properties of the Universe 
(Sect. |33t , and baryon statistics (Sect. [3~l4l in different models. 



3.1 Structure evolution 

We start by showing the basic features of the simulations. In Fig.Q~| 
we give a pictorial representation of the simulation with side of 
100 Mpc//i, and display the mass-weighted temperature evolution 
for a slice centered in the middle of the box, along the z-direction, 
having a thickness of ~ 7Mpc//i, and smoothed over a grid of 
1024x1024 pixels. We consider the /nl=0 (first two rows) and 
the /nl=1000 (last two rows) cases, and highlight the differences 
among them. At high redshift, z > 20, the /nl=1000 case shows 
more perturbations with respect to the /nl=0 one, and this is vis- 
ible at 2 ~ 10 — 15, as well, when feedback effects from early 
star formation heat the medium above ~ 10 4 — 10 5 K. The maps 
of the /nl=1000 case present more advanced stages of structure 
formation, and this keeps the gas slightly hotter, with temperatures 
of ~ 2 x 10 6 K, until z ~ 6, i.e. for about the first Gyr of the 
Universe. Temperature variations are up to a factor of < 2 for the 
/ NL =1000 case, but only up to ~ 10 per cent, for the other ones. 
Later on, at z < 6, more differences are not clearly visible and the 
final trends at lower redshift (z ~ 1) catch up and lead to the same 
behaviour at z ~ 0. 

We conclude that baryon history is affected by primordial non- 
Gaussianities mostly at early times, and this statement is strength- 
ened from Fig. [2] where we plot the molecular fraction at redshift 
z = 25, for the high-resolution runs, with side of 0.5 Mpc/h. We 
compare the /nl=0, the /nl=100, and the /nl=1000 cases during 
the very first stages of structure formation via molecular (mainly 
H2 and HD) cooling and catastrophyc run-away collapse. While 
in low-density environments there are no strong net differences, in 
the high-density regions the molecular content can change up to a 
factor of a few. Indeed, in the / NL =100 case, molecules increase 
of ~ 10 per cent (more exactly of 6 per cent) with respect to the 
/nl=0 one, and in the /nl=1000 case, of a factor of ~ 2.5. 
In Fig. [3] we highlight the differences by plotting the contrast 
between the / NL =100 and /nl= cases (top), and between the 
/ NL =1000 and /nl=0 cases (bottom), defined, for each pixel in the 
map, as 



C n 



^molecules, 1 ^molecules, 
^molecules,! H~ ^molecules, 



(2) 



In the previous equation, a; mo i e cuics,i is the molecular fraction in 
the pixel for the / NL =100 or the /nl=1000 case; while a; m oiecules,o 
is the reference /nl=0 case. It is clear that the comparison between 
/ NL =100 and /nl=0 expects much smaller differences than the 
one between /nl=1000 and f NL =0. Indeed, the top panel presents 
about 54 per cent of the entire area corresponding to contrasts of 
< 10 -3 (blue), almost 45 per cent corresponding to contrasts of 
~ 10 -3 — 10~ 2 (cyan-green), and the remaining few per cent reach- 
ing values of ~ 10 _1 , or larger (yellow-red). On the other side, the 
bottom panel displays larger deviations from the reference /nl=0 
case, with only 15 per cent of the region having contrasts < 10 -3 
(blue), roughly 70 per cent of ~ 10 -3 — 10~ 2 (cyan-green), and 
the remaining 15 per cent larger values (yellow-red). This means 
that larger /nl will induce higher molecular production, stronger 
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Figure 3. Molecular fraction contrast at redshift z = 25 between the 
/nl=100 and /nl=0 cases (top), and between the /nl=1000 and /nl=0 
cases (bottom). The color scales refer to the logarithm of the absolute value 
of the molecular fraction contrasts in the two cases. 



cooling, earlier gas collapse and earlier star formation, as we will 
see in the next section. 

The fundamental reason why gas and baryonic structures are sen- 
sitive to different /nl is that primordial non-Gaussianities affect 
the growth and evolution of the underlying dark matter haloes. As 
a consequence, the resulting baryon history reflects the imprints of 
the primordial dark-matter non-Gaussian distribution. 



3.2 Star formation 

At this point, we study the star formation history in the different 
models, paying attention to both the popIII and the popll regime. 
In Fig. [4] we plot the star formation rate densities from the sim- 
ulations of the 0.5 Mpc/h side boxes (left), and 100 Mpc//i side 
boxes (right), with the relative popIII contribution (lower panels). 
The simulations have similar trends, but the onset of star forma- 
tion is slightly earlier for the runs in 0.5 Mpc/fr side boxes. This 
is essentially due to the fact that the higher-resolution simulations 
are able to capture better molecular cooling and gas collapse in pri- 
mordial mini-haloes. Thus, they can address star formation in those 
regimes (i.e. at masses smaller than ~ 10 s Mq / h) where the larger 
100 Mpc//i side boxes are limited by resolution (see properties in 
Tab. Q}. We note that at both small and large scales, the main ef- 
fect of non-Gaussianities is a general shift of star formation, for 
the two population regimes. This is due to the fact that at larger 
/nl the contribution to the exponential tail of the mass function 
increases, leading to higher and earlier star formation. The differ- 
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Figure 1. Mass-weighted temperature evolution for the large-scale 100 Mpc//i side boxes, with /nl=0 (first and second row), and /nl=1000 (third and fourh 
row), at different redshifts (see legends). All the maps refer to slices centered at the mid-plane of the box, height z = 50Mpc//i, and have thickness of 
~ 7 Mpc / h, i.e. 1/14 the boxsize. 
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2 =25.00 



I =25.00 






—A ^ 
















f 1 








l-e.o 1 




-5.02 



/nl=0 



Atl=100 



Figure 2. Molecule maps at redshift 2 = 25 for the high-resolution 0.5 Mpc/ftside boxes, with /nl=0 (left), /nl=100 (center), and /nl=1000 (right). All 
the maps refer to slices centered at the mid-plane of the box, height z = 0.25 Mpc/ft, and have thickness of ~ 0.03 Mpc/h, i.e. 1/14 the boxsize. 



total 



L = 0.5 Mpc/h 




total 



L = 100 Mpc/h 



25 




20 



25 



Figure 4. Star formation rate densities as a function of the redshift for the simulation runs of the 0.5 Mpc/h side boxes (left) and 100 Mpc//i side boxes 
(right). In the top panels, the upper lines refer to the total star formation rate densities, while the bottom lines to the popIII ones. In all the cases, we consider 
/nl=0 (black solid lines), /nl=10 (cyan dotted lines), /nl=50 (green short-dashed lines), /nl=100 (red long-dashed lines), and /nl=1000 (blue dotted 
short-dashed lines. They correspond to Run05.0, Run05.10, Run05.50, Run05.100, Run05.1000, on the left panels, and to RunlOO.O, RunlOO.lO, Runl00.50, 
RunlOO. 100, RunlOO. 1000, on the right panels. Below, we plot the corresponding popIII contributions to the total star formation rates. 



ences among the various onsets are very small for /nl=0, /nl=10, 
/ NL =50 and /nl=100, and correspond to some ~ 10 7 yr. Only 
for the /nl=1000 case there is a larger gap of Az ~ a few, i.e. 
~ 5 x 10 7 yr, in the 0.5 Mpc/h box, and ~ 10 8 yr, in the 100 Mpc//7 
box. The behaviours of the popIII regimes are very little affected 
and their overall contributions drop down to ~ 10~ 3 — 10~ 2 , for 
any /nl, by redshift z ~ 10. Due to the anticipated evolution of the 
runs with larger /nl, the corresponding popIII contributions drop 
down earlier and stay slightly below the /nl=0 cases. On small 
scales (left panels), the difference is about ~ 2 orders of magni- 
tude, at z ~ 19 — 20, but only a factor of a few at later stages. 
At z ~ 10, in the /nl=0 case one expects a popIII contribution 
of ~ 0.0010, while in the / NL =1000 case, of ~ 0.0008 - 20 per 
cent smaller. Similarly, at the same redshift, on larger scales (righ 
panels), the differences are well within a factor of 2, with a popIII 



contribution of ~ 0.0097 for /nl=0, and ~ 0.0055 for / NL =1000. 
To check for degeneracies among different physical processes in- 
volved in the description of cosmic structure evolution, we have 
studied also how the scenario changes in two additional cases: 

• by adopting a common Salpeter-like IMF for the primordial 
popIII generation; 

• by comparin g with second-order correction s to the linear per- 
turbation theory jTseliakhovich & Hiratal 12010). which are sup- 
posed to generate primordial gas bulk mot ions, to suppress star for- 
mation on scales smaller than a few Mpc dMaio et alj2011 ), and t o 
induce delays in the reionization of the Universe (Dal ai et al J2010h . 

In Fig. [5] we plot the star formation rate densities for the 100 Mpc/h 
side boxes with primordial / NL =0, /nl=10, /nl=50, /nl=100, 
and /nl=1000. In these cases (see details in Tab. [TJ, we have 
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considered a Salpeter-like popIII IMF that is expected to cause a 
longer popIII-dominated epoch, because of t he loger lifetime s of 
the ~ 10s Mq stars dying as supernovae (see lMaio et al.lEoiol) . In 
fact, we find that the popIII regime dominates at early times and its 
contribution drops down to ~ 10 per cent only after ~4x 10 s yr. 
Furthermore, the global popIII contributions are higher in compar- 
ison with the top-heavy popIII IMF previously discussed, and the 
difference is more than one order of magnitude. At redshift z ~ 10, 
it is roughly a few times 10~ , which means that the primordial 
IMF can have impacts on the baryon and stellar evolution that are 
much stronger than the non-Gaussian effects. Therefore, if we want 
to give constraints on primordial non-Gaussianities it will be impor- 
tant first to properly constrain the popIII IMF. 
About the second check, we use only the small-box simulations, 
since they are more suitable to take into account primordial stream- 
ing motions. These are predicte d to have rms velocities of the order 
of ~ 30km/s at decoupling (Tseliakhovich & Hiratal l201ol) . So, 
we focus (see details in Tab. QJ on the same initial co nditions de- 
scribe d above for /nl=0 and /nl=100, and add (as in lMaio et al.l 
1201 lb a primordial bulk shift of 30 km/s to the gas. To consider sta- 
tistical deviations from the rms value, we additionally run the cases 
with 60 km/s and 90 km/s shifts (i.e. 2 and 3 times the rms value), 
as upper limits, which can be found in rare regions (less than ~ 1 
per cent) of the Universe, though. In Fig. [6] we compare the results 
for the two different /nl, and plot the /nl=0 and no-shift case 
(thin solid line) together with all the /nl=100 cases (thick lines). 
Obviously, star formation sets in at different times, according to 
the different values of the streaming motions and the resulting time 
delays are of the order of tens of Myr (Az ~ a few). However, 
if we compare the /nl=0 (thin solid line) and the / NL =100 (thick 
lines) cases for a fixed bulk shift velocity, we notice that the effects 
of non-Gaussianities are comparable to or negligible with respect 
to the impacts of the streaming motions. The star formation rate in 
the / NL =100 and 30 km/s-shift (thick dotted line) case almost coin- 
cides with the one in the /nl=0 run, and this leads to a degeneracy 
between these two different phenomena. Moreover, in the / NL =100 
and 60 km/s-shift (thick short-dashed line) case and in the / NIj =100 
and 90 km/s-shift (thick long-dashed line) case streaming motions 
cause much stronger delays. The value of 30 km/s is the average 
expected in the whole Universe, but different values can be found 
in different regions, thus, purely non-Gaussian effects are going to 
be strongly twisted with the statistical variations of such second- 
order contaminations. 

Finally, we stress that variations on the determination of ag of a 
few percents could produce a similar degeneracy, and different pri- 
mordial Lyman- Werner (LW) radiation strengths coming from the 
early star forming episodes could probably have even larger effects. 
In fact, the LW radiation (which is not explicitely followed in the 
simulations persented here) could dissociate primordial molecules, 
partially inhibiting H2 cooling in the neighbouring sites. Thus, the 
star formation process and the popIII formation at early times are 
expected to be delayed to lower redshifts. But the metals which are 
then ejected from primordial stars are insensitive to LW radiation: 
independently from molecule formation, they strongly enhance the 
cooling capabilit ies of the surrounding gas at any temperature (e.g. 
iMaio et alj2007b. and ca n rapidly lead to a popll-dominated regime 
dMaio et al. 2010. 1201 lb . Given the tight connections between pri- 
mordial IMF, LW radiation and chemical feedback more detailed 
studies will be needed to draw definitive conclusions on this issue, 
but definitely these will not change the results related to the /nl 
effects. 
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Figure 5. Star formation rate densities as a function of the redshift 
for the simulations of the 100 Mpc//i side boxes, with Salpeter popIII 
IMF. In the top panels, the upper lines refer to the total star for- 
mation rate densities, while the bottom lines to the popIII ones. In 
all the cases, we consider /nl=0 (black solid lines), /nl = 10 (cyan 
dotted lines), /nl=50 (green short-dashed lines), /nl=100 (red long- 
dashed lines), and /nl=1000 (blue dotted short-dashed lines. They cor- 
respond to RunlOO.O.SL, RunlOO.10.SL, Runl00.50.SL, RunlOO.100.SL, 
Run 100. 1000. SL. Below, we plot the con'esponding popIII contributions to 
the total star formation rates. 
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Figure 6. Upper panel: star formation rate densities as a function of the 
redshift for the simulation runs of the 0.5 Mpc/h side boxes with /nl=100. 
and primordial gas bulk velocities of km/s (thick black solid line), 30 km/s 
(thick cyan dotted lines), 60 km/s (thick green short-dashed lines), and 
90 km/s (thick red long-dashed lines). As a comparison, we overplot the re- 
sults for the 0.5 Mpc//i side box with /nl=0 and km/s for the primordial 
gas bulk velocity (thin black solid line). Bottom panel: the corresponding 
popIII contributions to the total star formation rates. 
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redshift 

Figure 7. Redshift evolution of the average pollution in the enriched regions 
(upper lines), and the average pollution in the whole box (lower lines), for 
/nl=0 (solid lines), /nl=100 (short-dashed lines), and /nl=1000 (long- 
dashed lines), in the 100 Mpc//i side boxes. The horizontal dot-dot-dot- 
dashed line is drawn at the critical value of 10 — 4 Z@. 

3.3 Chemical evolution 

We display, in Fig. [7] the metallicity evolutions for the runs of 
the /nl=0, /nl=100, and / NL =1000 cases, in the 100 Mpc//j side 
box with a top-heavy popIII IMF. The upper lines correspond to 
the average enrichment (i.e. the average metallicity of the pol- 
luted regions only), and the lower lines to the average metallic- 
ity (in the whole simulated volume), as a function of the redshift. 
What is particularly e vident is the rap idity of early metal pollu- 
tion (consistently with lMaio etalfc o 10). Indeed, the critical value, 
Zcrit = 10 -4 ZQ, is reached in only a few 10 7 yr, from the first 
episodes of star formation and the primordial yields are so high 
that effi ciently pollute the surrounding medium (see more discus- 
sion in iMaio et alj|20ich . The following metal evolution shows 
an increasing trend that leads to an average enrichment of about 
Z ~ 10~ 3 - 10~ 2 , i.e. ~ 0.1 Z Q . Due to the patchyness of the 
metal enrichment process, it is possible to get very strong local 
pollution, too, with highly supersolar Z. The average metallicity in 
the whole volume (solid line) increases more gradually, by chang- 
ing of ~ 6 orders of magnitude between z ~ 16 and z ~ 6. This 
behaviour is due to the fact that, as cosmic time evolves, the re- 
gions involved in metal pollution are larger and larger, in compar- 
ison to the very first isolated events. However, the average enrich- 
ment reaches the critical value, Zcrit, by redshift z ~ 6, when 
the Universe is roughly one-billion-year old. Afterwards, during 
the subsequent ~ 13 Gyr, Z increases of ~ 3 orders of magni- 
tude, and the cosmic star formation history is fully dominated by 
the popll regime (see also Fig. |4j. Among the different elements, 
oxygen is the dominant one, since it is heavily produced by PISN 
and SNII explosions. While at early times it dominates by about 
one order of magnitude, at later stages other elements become im- 
portant, in particular, iron, which is mostly produced by the death 
of long-living stars. It catches up carbon and overtakes a-elements, 
like magnesium, sulphur, and silicon. 

Briefly speaking, the fact that the average enrichment reaches the 
critical value already by z ~ 6 has broad implications also for 



reionization, given that metallicity affects the spectral energ y dis- 
tributions (SEDs) of stellar populations (e.g. ISchaererll2003h . In 
particular, primordial stellar populations have UV fluxes (at wave- 
lengths shorter than 912 A) which are up to four times larger than 
the corresponding popll ones, and in principle, could ionize hydro- 
gen in a much more efficient way. However, given the high level 
of enrichment in the clustered star forming regions (upper lines in 
Fig. [7} the popIII star formation rate at z ~ 6 is only roughly 
~ 10" 2 - 10" 4 the total one (depending on the modeling). The 
relative number fraction of popIII hundred-solar-masses stars with 
respect to the corresponding popll stars drops further down of a 
factor of ~ 10 — 100. As a consequence, it is unlikely that popIII 
stars will dominate the reionization process simply because they 
are too rare and, additionally, shine for a too short time (up to few 
10 6 yr) to provide significant amounts of UV photons. 
T hese conclus i ons co nfirm and extend to lower redshift the finding 
bv lMaio et alj J20101) . and show how primordial non-Gaussianities 
affect the whole picture. In fact, by comparing the different /nl 
cases, it emerges that the strongest difference is found for the 
/ NL =1000 run, in which metal enrichment kicks in at earlier times, 
tracing back the star formation rate (see Fig. |4]. Its behaviour is 
not drammatically different from the other ones, below z < 12, 
though. We note that the /nl=0 and / NL =100 runs predict ex- 
tremely similar evolutions, with average metallicities differing by a 
factor of a few only at the very beginning of the pollution process 
(at z ~ 13 — 15, where the average enrichment is of Z < 10 10 ), 
and rapidly converging later on. 

The trends found in all the other runs show very similar behaviour 
and lead to the same conclusions. 

3.4 Baryon distributions 

After discussing the general behaviours of the baryons in the Uni- 
verse, it is interesting to go through their clumping properties and 
to investigate how the resulting "visible" objects are affected. So, 
we identify gas clouds and dark-matter haloes by means of a FoF 
algorithm, as already described in Sect. [2] 

To probe the effects of non-Gaussianities on the small primor- 
dial gas clouds, in Fig. [8] we show the cumulative distributions at 
early times for the 0.5 Mpc//i side boxes. We focus on the /nl=0, 
/nl=100, and /nl=1000 cases, at redshift z = 25 and z — 15. 
At z ~ 25, the primordial non-Gaussian perturbations have well 
visible consequences in the gas dumpiness and lead to more mas- 
sive and more numerous clouds with respect to the standard Gaus- 
sian scenario, for gaseous structures with mass < 10 7 Mq//i . The 
differences are about a factor of ~ 2 — 3 in the /nl=1000 case, 
but only ~ 10 per cent for the / NL =100 case. This reflects the 
larger probability at higher /nl of forming objects that are able 
to trap and condense primordial cosmic gas. The following growth 
of structures washes out the differences in a few ~ 10 8 Myr, and 
by z ~ 15 the gas distributions for these small objects are almost 
undistinguishable. 

The following evolution will lead to the formation of protostel- 
lar cores, to the birth of the first generation of stars, and then to 
the establishment of the standard stellar populations (see previous 
discussion). To account for the differences which persist down to 
lower redshift on the stellar evolution properties, we plot in Fig. [9] 
the stellar mass distributions in the range ~ 10 7 - 10 13 M Q , at 
redshift z = 8.51 and z = 3. We consider the / NL =0, /nl=100, 
and /nl=1000 cases in the large boxes, because they have higher- 
mass samples. We find that also the stellar distributions are af- 
fected, mostly at earlier times, with differences that can reach a 
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Figure 8. Cumulative primordial-cloud distributions as a function of gas 
mass, for z = 25 (upper panel) and z = 15 (lower panel), for the runs 
with /nl=0 (black solid lines), /nl=100 (red long-dashed lines), and 
/nl=1000 (blue dotted short-dashed lines). We plot results for the high- 
resolution simulations with 0.5 Mpc/ft side boxes, labelled as Run05.0, 
Run05.100, Run05.1000, respectively. 



Table 2. Observational determinations of the /nl parameter, according to 
different authors and different techniques. 



Reference 


Range for / NL 


Confidence 


Komatsu et al. (2002) 


[-3500, +2000] 


2(7 


Komatsu et al. (2003) 


[-58, +134] 


2(7 


Gaztanaga & Wags (2003) 


[+4, +18] 


1(7 


Gaztanaga & Wagg (2003) 


[-4, +5] 


1(7 


Chen & Szapudi (2005) 


[-260, +40] 


1(7 


Chen & Szapudi (2006) 


[-30, +74] 


1(7 


Cabella et al. (2006) 


[-80, +80] 


1<T 


Cabella et al. (2006) 


[-160, +160] 


2(7 


Spergel et al. (2007) 


[-54, +114] 


2(7 


Afshordi & Tolley (2008) 


[+109, +363] 


2(T 


Slosar et al. (2008) 


[-29, +70] 


2(T 


Slosar et al. (2008) 


[-65, +93] 


3(7 


Slosar et al. (2008) 


[-31, +70] 


2a 


Slosar etal. (2008) 


[-96, +96] 


3(7 


Hikage et al. (2008) 


[—70, +91] 


3(7 


Yadav & Wandelt (2008) 


[+27, +147] 


2(7 


Komatsu et al. (2009) 


[-9, +111] 


2(7 


Smith et al. (2009) 


[-4, +80] 


2<T 


Curto et al. (2009) 


[-8, +111] 


2a 


Curto et al. (2009) 


[-18, +80] 


2a 


Vielva & Sanz (2009) 


[-94, +154] 


2a 


Rudiord et al. (2009) 


[4-44, +124] 


2a 


Rudiord et al. (2010) 


[+42, +104] 


2a 


Hou et al. (2010) 


[+34, +120] 


1(7 


Hou et al. (2010) 


[+12, +82] 


la 


Cabella et al. (2010) 


[-9, +85] 


la 


Komatsu (2010) 


[+11, +53] 


la 


Komatsu (2010) 


[-10, +74] 


2a 




7 8 9 10 11 12 13 14 
Log 10 ( M 5lor [M /h] ) 



L = 100 Mpc/h 


'nl = 

'nl = 100 

f NL = 1000 


z = 3.00 


L-i L. 

n , i 



7 8 9 10 11 12 13 14 
Log 10 ( M slor [M e /h] ) 



Figure 9. Cumulative cloud distributions as a function of stellar mass, 
for z = 8.51 (upper panel) and z = 3.00 (lower panel), for the runs 
with /nl=0 (black solid lines), /nl=100 (red long-dashed lines), and 
/nl=1000 (blue dotted short-dashed lines). We plot results for the high- 
resolution simulations with 100 Mpc//i side boxes, labelled as RunlOO.O, 
RunlOO.100, RunlOO.1000, respectively. 



factor of ~ 10 for the /nl=1000 case and a factor of a few for the 
/nl=100 case, in the high-mass tail (> 10 9 M s IK). At the low 
mass-end (~ 10 7 — 10 s Mq / h), instead, differences are only a 
factor of ~ 2 for /nl=1000, and ~ 10 per cent for / N l=100. At 
low redshift, z ~ 3, some differences persist for masses larger than 
~ 10 11 Mq/Ii, but the most of the objects shows well converg- 
ing trends. Below redshift ~ 3, the residual differences are further 
washed out by the subsequent structure formation and the on-going 
mass assembly process, so at z ~ — 1 there are no distinctions at 
all in the various stellar components. 



4 DISCUSSION AND CONCLUSIONS 

Cosmological structure formation is strongly affected by the ini- 
tial perturbations which it originates from. As a consequence of the 
central-limit theorem, it is common to assume a Gaussian distri- 
bution, but detailed studies of the CMB leave some room for de- 
viations from the expected behaviour, and for an /nl parameter 
different from zero. Since non-Gaussianities alter the "initial con- 
ditions" of the problem of structure formation, and, thus, influence 
the temporal evolution of the forming structures, we have addressed 
their impacts on the visible, baryonic matter. 
Given the lack of extensive works on this subject, we have per- 
formed the first high-resolution cosmological N-body, hydrody- 
namical, chemistry simulation of cosmic structure formation, in 
the frame of the ACDM model and in presence of primordial 
non-Gaussianities. We have studied the principal properties of 
gas and stars in different scenarios, by following: cosmic evo- 
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lution of e~ H, H+, H~ He He+, He++ H 2 , H+ D, D 
HP, HeH+ dYoshida et all 120031 ; Im^o et alj 120061. I2001J 200 1 
l2010h : radiative gas cooling at high dSutherland & Dopitalll99: 
and low temperatures (both f rom molecules a nd fine-structure 
metal tran sitions, according to Maio et al.ll2007t) : multiphase star 
formation dSpringel & Hernquistll2003l) : UV background radiation 
dHaardt & Mada uj 1 19961); wind feedback dSpringel & Hernquistl 
120031 : lAguirre et all l200lh : metal pollution from popIII and/or 
popll stellar generations, ruled by a critical metallicity thresh- 
old of Z C rit = 10" 4 Z dTornatore et aill2007t iMaio et al.|[2oToh 
and leading gravi tational enrichment into the surrounding medium 
dMaio etalj|2010h . 

Observational evidences and continuous investigations over the last 
decade (see details in Tab. [2]l have progressively constrained /nl 
to positive values, thus, we have considered /n L =0, 10, 50, 100, 
and 1000. From Tab. [2] the most likely valued are expected to be 
between /nl~ — 100, and the /nl=1000 case has to be simply 
considered as an extreme upper limit. 

The goal of this work has been to theoretically explore the effects 
of different /nl values on the structure formation properties and 
on the baryon history. Thanks to the high resolution reached (up to 
~ 40 Mq /ft), we are able to address the cosmic evolution from the 
very early times - when the first star formation episodes take place 
- till to the present. 

Our main results (Sect. [3} have shown that structure formation pro- 
cesses (Sect l3. lb . star formation history (Sect. 13.21 . chemical evo- 
lution (Sect. [3~3l >, and visible-matter distributions (Sect. [3~l4l are all 
affected by non-Gaussianities. However, while the effects are quite 
evident for the extreme /nl=1000 case, they are much milder for 
smaller, more ralistic values of /nl. 

Indeed, gas temperature and composition reveal some effects of pri- 
mordial non-Gaussianities only at very high redshift, with changes 
of the order of ~ 10 per cent, at z > 10, for /nl=0 — 100 and of a 
factor of ~ 2 for / NL =1000 (Fig.Q)- 

Primordial molecular fractions (Fig. |2]and Fig. [3} are particularly 
senstive to the thermodynamical properties of the medium and their 
behaviour increases mostly in the high-density regions, where the 
effects on clustering due to primordial non-Gaussianities are more 
relevant. 

Star formation history (Fig. |4} has a converging behaviour despite 
of non-Gaussianities. The main discrepancies among the models 
are found for the onset, which takes place at earlier epochs for 
larger /nl- The time-lags are around some 10 7 yr, at z ~ 15, and 
they are reflected in the transition from the popIII regime, described 
by a top-heavy IMF, to the popll regime, described by a Salpeter 
IMF (Fig.|4j. Despite such discrepancies, the popIII contributions 
stay around ~ 10~ 3 — 10~ 2 at z ~ 10 for all the models, with only 
~ 1 — 10 per cent scatter. 

Given the uncertainties on the primordial popIII IMF we also run 
corresponding simulations with Salpeter-like popIII IMF and found 
that (Fig. [5} the resulting deviations of the popIII star formation 
rates could be larger than the non-Gaussian ones. 
Additionally, we compare d with the second-order corr ections in 
linear perturbation theory dTseliakhovich & Hirat a 2010) and their 
impact on structure formation evolution dMaio et alj|201lb . show- 
ing that (Fig. |6) also the resulting primordial gas bulk flows 
could have comparable or larger repercursions on the visible mat- 



3 In Tab.ff] the CMB notation is assumed, thus, to compare with the values 
adopted in the simulations (/nl=0> 10, 50, 100, 1000), one has to multiply 
the tabulated data by ~ 1.28. 



ter. Therefore, the behaviour of baryons seems tightly bound and 
twisted with several physical mechanisms, whose effects are de- 
generate with non-Gaussianities. 

As a consequence of the trends of the star formation rates, metal 
enrichment results (Fig. |7J are affected by the initial /nl, mostly 
at early times, with pollution histories converging at low redshift 
(z < 12) and overcoming, in average, the critical metallicity at 
z ~ 6, when the Universe has an age of roughly 1 Gyr. This is quite 
independent from the primordial non-Gaussianities, popIII IMF, or 
additional second-order effects, but at high redshift metal spreading 
is highly sensitive to all these phenomena and can show degenera- 
cies with them. 

Gaseous and stellar mass statistics (Fig. [8] and Fig. [9} show vari- 
ations at the high-mass tail, and rare, more massive objects are 
more common in higher- /nl cosmologies. The distributions of the 
very first structures show an increment of a factor of a few in the 
/nl=1000 case, but only some ~ 10 per cent for /nl= — 100, 
mainly at the high-mass end. However, the subsequent evolution at 
later times washes them out almost completely. 
Throughout this work, we have used the /nl parameter only to 
quantify deviations from the Gaussian shape. We have not inves- 
tigated models with a negative /nl as they have been progres- 
sively dismissed by observations. However, for sake of complet- 
ness, we note that, theoretically, one would expect roughly specu- 
lar behaviours with respect to the /nl=0 case when universes with 
negative /nl are considered - a result already well established from 
simple dark-matter-only simulations (as mentioned in Sect.QJ. In- 
deed, in such scenarios, the probability distribution function of the 
initial density field would be skewned towards the underdense tail, 
causing a paucity of the more massive objects and, in general, a 
delayed evolution of the cosmic structures. We therefore expect the 
variations from the /nl=0 case originated from a positive /nl pa- 
rameter to be accordingly reversed when the corresponding nega- 
tive /nl parameter is adopted. 

Besides the /nl formulation, higher-order parameters can be in- 
troduced, but their unc ertainties are hug e with, e.g., —7.4 x 
10 5 <3nl< 8.2 x 10 s dSmidt et aljboiol) . Thus, a study of their 
impacts on the cosmic structures would be quite vague and elu- 
sive. Furthermore, the effects of the second-order /nl parameter 
are quite small, so it is reasonable to think that inclusion of higher- 
order parameters will not dramaticaly change our conclusions. 
Even if future initial conditions will include additional higher-order 
corrections, being definitely more complete and detailed, probably 
they would not alter the global scenario we have drawn. For what 
we said in the previous sections, other physical phenomena, like 
cosmic rayfl radiative transfer and reionizatior0, magnetic field^fl 



4 Se e e.g. Ginzburg & Svrovatski il dl964l); |Mannheim & Schl ickeiseil 
<1994l): iBiermann & Sigll {200lh; ISigl^Tatl 120041); IStane-vl | 20o X 

l) ;|jascheet al.l J2007I) 



Vasiliev & Shchekinov 1 2006 l); |Pfrommer et al. Id2007l);|jasche et al.H200' 
5 See e.g . [Abel et alj jl999l): iGnedin & AbeT feOOll); ICiardi et al 
|200lL 120031); llVlaselli et alj 120031); IWhitehouse et alj d2005t); llliev et al 
{2006 1 ) ; [Mellema et alj l2006h ; llrac & Cenl 120071); lAubert & Tevssiei 

12008) ; jjohnson et alj <2008t); iPawlik & Schavei j2008h; iFinlator etaT 

12009) ; IPetkova & Springell 12009. 2010); ICantalupo & Porcianil d201ll) 



IWise & Abell feOlCl). 



B See e .g lOrszag & Tand dl97Sj); 
Brio&Wj dl988l); bolag et alj ll999l); 



Londrillo & Del Zannal l2000l); lT6thT 



Phillips & Monaghai 



Briig gen et alj l2005l); iPrice & Monaghanl ilOoSi 
Stone et alj <2008|); livu et alj 120081); Ipolag & 



Brandenburg (2010); 
KotarbaetallfeOld) . 



Balsara & Spice; 
20001); lB0rveetal 
Li et al. 
Stasyszyn 

Biirzle et al.l d2010l) ; [Gaburov & Nitadori 



1985) 



1999); 



2001); 
42008); 
12009 ); 
(2010); 



Baryon history and non-Gaussianities 1 1 



AGN feedback etc., are not supposed to experience strong modi- 
fications, because of primordial non-Gaussianities. 
We stress that even small variations in the determinations of the 
cosmological parameters, above all errors of a few per cent in 
the estimates of as, could add further degeneracies to disentangle 
the purely non-Gaussian results from the primordial IMF, gas bulk 
flows, etc. 

We summarize our conclusions by stating that for reasonable values 
of /nl= 0—100 in the primordial non-Gaussianities (in agreement 
with observations): 

• early molecular fractions in the cold, dense gas phase are pre- 
dicted to be altered of ~ 10 per cent; 

• small temperature fluctuations of < 10 per cent are induced 
during the first Gyr of the cosmic evolution; 

• the onset of the first star formation events, at z > 15, is influ- 
enced of some 10 7 yr; 

• the epoch of poplll/popll transition is affected by some 
10 7 yr; 

• variations of < 10 per cent are found in the mass distributions 
during the formation of first structures; 

• only mild variations in the chemical history of the Universe 
are expected; 

• extremly high values for /nl would imply corresponding vari- 
ations of a factor ~ 2 — 3 in the thermodynamical and chem- 
ical properties of the gas, earlier star formation histories of ~ 
0.5 — 1 x 10 8 yr, and changes up to a factor of ~ 10 in the high- 
mass end of the gas cloud and stellar mass distributions at early 
times; 

• non-Gaussian effects could have degeneracies with other 
physical processes, like primordial gas bulk flows, with unknown 
star formation quantities, like popIII IMF, or with the determination 
of cosmological parameters, mostly as at a few per cent level. 
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